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Abstract 



We consider AdS2 x R 2 solutions supported by a magnetic field, such as those which arise in 
the near-horizon limit of magnetically charged AdS^ Reissner-Nordstrom black branes. In 
the presence of an electrically charged scalar field, such magnetic solutions can be unstable 
to spontaneous formation of a vortex lattice. We solve the coupled partial differential equa- 
tions which govern the charged scalar, gauge field, and metric degrees of freedom to lowest 
non-trivial order in an expansion around the critical point, and discuss the corrections to 
the free energy and thermodynamic functions arising from the formation of the lattice. 
We describe how such solutions can also be interpreted, via S-duality, as characterizing 
infrared crystalline phases of conformal field theories doped by a chemical potential, but in 
zero magnetic field; the doped conformal field theories are dual to geometries that exhibit 
dynamical scaling and hyperscaling violation. 



Contents 



1 Introduction [3] 

2 Magnetic AdS 2 X i? 2 solutions H 

3 The vortex lattice U 

3.1 Basic droplet [7] 

3.2 Vortex lattice E 

3.3 The radial equation and boundary conditions [TO] 

3.4 Higher order corrections to the gauge field and metric [12] 

3.5 0(e 2 ) stress-energy at the wall H5] 

3.6 Pictures of the modulated phase [16] 

4 Comments on physics of the lattice model 1181 

5 Connecting with more general gravity solutions 1191 

5.1 Basic EMD theory and magnetic solutions [20] 

6 Discussion [22] 



2 



1 Introduction 



A topic of recent interest has been the holographic description of phases of quantum field 
theory with spatial anisotropy and/or inhomogeneity (T 15 . This is motivated in part 



by the crucial role that momentum relaxation due to inhomogeneities plays in transport 
phenomena in condensed matter systems, and in part by intrinsic interest in the rich physics 
of such phases. 

Our goal in this work is twofold. On the one hand, as an extension of the ideas discussed 
in 16 , we would like to illustrate the emergence of crystalline ground states ('solids') 
in conformal field theories doped by a chemical potential coupling to a globally conserved 
U(l) charge, but in zero magnetic field. In 2+1 dimensions, monopole operators associated 



with the global U(l) symmetry 16 - 18 serve as order parameters for solid phases in doped 
CFTs. Electric-magnetic duality allows one to find a dual description where the magnetic 
degrees of freedom are manifested in terms of electrically charged operators. In the bulk 
gravitational description, this allows us to view the formation of the solid by studying vortex 
lattice formation in the theory of a charged scalar moving in a background magnetic field. 
A crucial advantage of studying the solid phases of doped CFTs by using this dual charged 
scalar is that the Dirac quantization condition on the monopole charge translates into an 
exact commensurability relation between the area of the unit cell of the crystal and the 



density of doped charges 16 19] 



On the other hand, an open problem in the study of holographic lattices has been to find, 
analytically, gravitationally back-reacted solutions for a crystalline lattice of dimension 
d > 1 . This has largely been because of the relative difficulty of solving coupled systems of 
partial differential equations, instead of the ordinary differential equations which normally 
govern simple backgrounds in gauge/gravity duality. Here, we give an example of such 
a crystalline metric in d = 2. Our work builds on the earlier papers [3], which found 



an elegant solution for a vortex lattice in the probe approximation, and 11 , where the 



backreaction of such a lattice on bulk gauge fields was studied in a different setting. 

The organization of this paper is as follows. In §2, we review the basic unperturbed 
AdS 2 x R 2 solution. In §3, we incorporate a charged scalar field and describe the vortex 
lattice solution. In §4, we describe the basic physics visible in the perturbative vortex 
lattice solution. In §5, we characterize how such a lattice could also emerge in the IR 
geometry of a gravitational solution which exhibits dynamical scaling with hyperscaling 
violation, as the S-dual of a doped CFT in zero magnetic field. Possible directions for 
future research are discussed in 56. 
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2 Magnetic AdS2 X i? 2 solutions 



Consider the theory with action 

S = 

where A < 0. It has AdS2 x R 2 solutions with metric 



r/.s 2 -= L 2 ( + + cfe 2 + rfy 2 



where 



1 

L2 



-2A . 



The gauge field supporting the solution is 

F xy = Q m dx A dy 
with Q m fixed in terms of the AdS radius by the equation 

Q m = V2L . 



(2.1) 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



In particular, this means that for these solutions, fixing the magnetic field fixes also the 
cosmological constant and the AdS radius. 

In addition to its intrinsic interest, this solution arises as the near- horizon geometry of 
extremal magnetically charged AdS/Reissner- Nordstrom black branes with AdS^ asymp- 
totics. In this context, the AdS% near-horizon region has played a crucial role in elucidating 



the non- Fermi liquid behavior of probe fermions 20 22 scattering off the bath of locally 



critical excitations represented by the AdS2 geometry 23 24 . 



3 The vortex lattice 



Our interest is not in the pure AdSz solution (2.2). We wish to include also an electrically 
charged scalar field, ip, in the full action. In part, this is because a generic such theory 
could include such scalars; in part, it is motivated by the duality considerations to be 
described in §5. 

In any case, here, we will see that in some ranges of parameters, the charged scalar will 
qualitatively change the IR physics. The simplest case in which we can see this effect will 
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be directly in the AdS 2 x R 2 background of § 2. We will impose a hard wall cutoff at r = r 
in the deep IR, along with suitable boundary conditions, to be described below. We can 
think of r as a proxy for a 'confinement scale' or a 'temperature.' Tuning the magnetic 
field relative to the 'temperature' will trigger the scalar instability. 

After including the ip coupling to the gauge field the action becomes 

S = J d'x^g (r-2A- \f 2 - \V^\ 2 - ™ 2 H 2 - A|Vf) • (3.1) 

We have defined V M = + ieA^, where e is the electric charge of the scalar field. From 
this point on we will set e = 1. This action has a stress-energy tensor of the form 

1 1 

T, u = -^£ mat + -F MA F* + e 2 ^M (3.2) 



where 
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C mat = -F 2 + |V^| 2 + mVl 2 + AH 4 - (3.3) 



We may expand the magnitude of IV^VI 2 as 

| V^\ 2 = \d^\ 2 + iA^r - ^9^) + A 2 H 2 . (3.4) 

At this point, we can calculate the Euler-Lagrange equation for ip by differentiating with 
respect to ip*: 

^(v^V^V) = -V^9(^V^ - m 2 i) - 2A|^| 2 ^), (3.5) 
and the equation of motion for the gauge field, 

! d ll {^-gF^)=i{i,d u r -rd v 4>) + 2A v \i,\ 2 . (3.6) 



'-9 



In addition to these equations of motion, we will also need to solve the Einstein equations, 

(.R-2A) 

when we include backreaction of the ip condensate on the gauge field and the metric. 

We will expand perturbatively in a small parameter e around the solution ip — with 
background gauge field of the form 

A x = Q c y,A y = 0. (3.8) 

The scalar field in the competing vortex phase will itself be of order e. For fixed r and 
boundary conditions (to be discussed below), we will choose Q c to be just at the onset for 
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the transition to forming vortices. At this critical value of the magnetic field, the ijj = 
solution will be degenerate with a vortex lattice solution. As we increase the magnetic 
field to slightly above its critical value, ip = will no longer be the preferred solution, and 
the vortex lattice will be preferred. As is familiar, the onset of the transition is signalled 
by the existence of a purely normalizable solution for ip that respects the IR boundary 
conditions. 

We can parametrize the backreaction of the scalar on the gauge sector through a pertur- 
bative expansion in the distance away from the critical field. The scalar will have the 
form 

V>(r, x, y) = e^i(r, x, y) + e 3 i/j 3 (r, x,y) + ... (3.9) 
and the gauge field will have the form 

A x (r,x,y) = Qy + e 2 a x 2 (r,x,y) + ... J (3.10) 
A y (r,x,y) = e 2 al(r, x, y) + . . . , (3.11) 

with A t = A r = 0. When we consider lattice solutions which are periodic in x and y, the 
backreaction of ip on the gauge field will require both A x and A y to be nonzero at 0(e 2 ), 
with both x and y dependence. 

A similar statement holds for the metric at 0(e 2 ). Our metric ansatz, to 0(e 2 ), will be 

ds 2 = L 2 j-^((— 1 + e 2 a(r,x,y))dt 2 + (1 + e 2 a(r,x,y))dr 2 ) + (1 + e 2 b(r, x, y))(dx 2 + dy 2 ) 

(3.12) 

Because at zeroth order in epsilon the AdS2 x R 2 metric is exactly supported by the 
magnetic field (i.e. the gauge field is not a probe), we find it necessary to include metric 
backreaction once we backreact on the gauge field. This distinguishes our situation from 



that considered in e.g. 11 



The radial magnetic field will be 

B r = Q + e 2 (d y a x 2 (r, x, y) - d x a y 2 (r, x, y)) (3.13) 



In general, when we backreact on the magnetic field, we may expect there to be a non- 
normalizable piece at order e 2 , i.e. A x (r — > 0) = (Q + 5Qe 2 )y. This shifts the naive critical 
value of the field at the transition. However, because the critical point is actually only 
dependent on the dimensionless combination Q/r 2 ., we can (and will) impose that there is 
no non-normalizable correction to the gauge field in our backreacted solutions. That is, 
we will set 5Q = 0. The value of the critical point will still have an 0(e 2 ) shift; it will 
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manifest itself as an 0(e 2 ) shift in the location of the hard wall, r$ — > tq + 5ro€ 2 . These 
two scenarios are equivalent; in both cases we should think of the backreaction of ip on the 
metric and gauge field as inducing a shift in the dimensionless parameter which controls 
the critical point at 0(e 2 ). 



3.1 Basic droplet 

We now examine the solutions of the field equation for ip, in the limit where we can neglect 
the back- reaction of ip on the gauge field and on the metric. (This will be at order e.) 
Very similar equations have been examined in the literature on vortices in holographic 
superconductors (3 25 26 . The basic building block for the solutions we will study is the 



"droplet" solution of 25 



We will begin by setting A = in the potential for the scalar and proceed with the metric 
and the mass of the (dualized) monopole field unspecified. Both of these will affect the 
radial solution for the scalar, but we will see that the spatial part of ipi decouples from the 
radial equation for all metrics we might consider, and so we can find the basic form of the 
droplet solution while leaving the metric general. 

For metrics with components which only depend on r and for which g xx = g yy , we can 
solve this equation by separation of variables, assuming that 

= PoP (r)g(y)e ikx , (3.14) 



where po is an overall constant. Inserting our choice of gauge (3.8) yields, after some 
algebraic manipulation, 

rr (^P» + -r^§z(V=99 rr )Pn(r)] - ^ (3.15) 
p n (r) \g xx y/—gg xx or J g xx 

= — t-^v {g'n(y) - (Qv + k ) 2 9n{y)) = -K, 
9n{y) 

where X n is the eigenvalue from the separation of variables. First we will consider the 
equation for g(y), which will yield the basic droplet solution. This solution will only exist 
in the parameter ranges which admit a normalizable solution to the radial equation; we 
will discuss this in the next section. The equation for g becomes, 

9'n-(Qy + k) 2 g n = X n g n . (3.16) 
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Now, redefining Y = y/Q(y + q), the g n equation becomes 



g'niX) - ( Y 2 + ^ ) 9n(Y) = 0. (3.17) 



Solving, we get that 



g n (Y) = c + D u+ (V2Y) + c_D v _ (iV2Y) (3.18) 
where c± are constants, v± — \ ^—1 ± and D v {x) is the parabolic cylinder function. 



The reader may recognize the differential equation for g n , (3.16), as the same eigenvalue 
problem that arises in the study of the quantum mechanics of the simple harmonic oscilla- 
tor. More properly, this is the case for appropriate choices of the separation constant. In 
these cases, we can write the (normalizable) solution for g n in terms of the familiar Hermite 
polynomials: 

g n = e~ Y2/2 H n {Y), (3.19) 

with eigenvalues X n = 2Q(n + 1/2). The nth eigenvalue here characterizes the nth Landau 
level of the ip particles. The "droplet" solutions with this shape were first discussed in the 



series of papers 25 , in a related but distinct context. The single droplet solution is when 
n — 0, which is just a Gaussian centered at y — —k/Q, g(y) = e~ y2//2 . Note that g n = 
constant is not a solution to the equations of motion. 



3.2 Vortex lattice 

Of more interest to us is a solution which preserves some discrete subgroup of the trans- 
lation invariance of the original system. The basic droplet of §3.1 breaks translations 
entirely. However, more symmetric solutions can be obtained by taking linear combina- 
tions of droplets, which still solve the (linearized) equations of motion neglecting back 
reaction. 

A vortex lattice can be constructed as follows |3j, using the zeroth Landau level solutions 
for the ip field. The basic solution is 

^ ( y; k) = e-£ = e-^ v+ i )2 . (3.20) 

An appropriate superposition to give a lattice in the x — y plane is 

1 oo 

*iat(z,y) = i c i elklX Mv, h) (3.21) 

l=— oo 
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where 



C = e , fe,= \[Q 3.22 



for arbitrary V\ and t>2- 

One can write this in terms of the elliptic theta function 6*3: 

00 

9 3 (v,r) = (1 = ^\ z = e iwv (3.23) 



{=— 00 

as 



ip!(x,y,r) = pop(r)Viat(x,y), ^i at (x,y) = e 2 9 3 (v,t) (3.24) 

with 



y/Q{x + ly) __ 2m - v 2 

, r = 

Vi vf 



(3.25) 



That the solution (3.21) represents a lattice is now evident from the basic properties of the 



elliptic theta function. For instance 

9 3 (v + 1,t) = 9 3 (v,t) (3.26) 

and 

e 3 ( v + t,t)= e- 2ni{v+T/2) 6 3 {v, r) (3.27) 

implying that lattice returns to its value (up to a phase) upon translation by the lattice 
generators 

a = -^uiSL b = -= (—d v + —d x A . (3.28) 

These have been fixed such that the area of a unit cell is 2tt/Q, containing exactly one 
flux quantum. It is this quantization condition which translates, in the electromagnetic 
dual, to the commensurability condition between the area of the unit cell and the density 



of doped charges 16 , 19 



That ^/i a t should be called a vortex lattice, despite the fact that it is composed of an array 
of the droplet solutions of [25], is evident from the fact that #3 vanishes on the lattice 
spanned by half-integral multiples of the lattice generators (giving rise to vortex cores), 
and has a phase rotation of 2n around each such zero. 

Some common lattice shapes are obtained by choosing particular values of the parameters 
vi,V2- A rectangular lattice can be obtained by setting t> 2 = 0. In this case all coefficients 



in equation (3.21) are equal, q = c = 1. The ratio of length to width of the rectangle 



is parametrized by v\. For the special choice V\ = V2tt, the lattice is square. Another 



9 



special choice is v% = \v\\ this yields a rhombic lattice. In this case q = 1 for I = mod 2 
and q = — i for I = 1 mod 2. For the special case ui = 2y/n the rhombus is square (but 
now rotated 45° w.r.t the x axis), and for V\ = the lattice is composed of equilateral 

3 3 

triangles (though the unit cell is still a rhombus). 

At this point, nothing has fixed the "moduli" t>i,t>2 of the vortex lattice, nor the overall 
magnitude po of if>\. In standard Landau- Ginzburg theories, apparently the triangular 
lattice is preferred. One could find preferred shapes in the approach here by including 
leading non-linearities in \if)\ in the free energy, and minimizing the free energy density. It 
might be interesting to do this, while introducing parameters to vary that could lead to 
phase transitions in the preferred lattice shape. 



3.3 The radial equation and boundary conditions 

Now we consider the differential equation for p(r). At order e we are still in an AdS2 x R 2 
background, which means that if) should scale as a power law in r. Choosing the solution 
that vanishes at the boundary, we get 

p(r) = r a (3.29) 

where a = § (l + ^1 + 4(Q + m 2 L 2 )) . 

At the hard wall cutoff, r = r , we will need to impose a consistent set of boundary 
conditions. One way to do this is to consider a method very similar to the prescription 



of 27 . We will add a mirror image of the spacetime to the other side of the wall and glue 
them together at the IR boundary, r = r . Thus, we have two asymptotic UV boundaries 
(in our coordinates at r = and r = 2r ) and mirror solutions for the metric and the 
fields on either side of the wall. We will require the metric and fields to be continuous at 
the wall, but their derivatives will have a discontinuity. That is, we impose Israel junction 
conditions at the wall, including any localized energy-momentum sources present there. 
At the end of the day, we can quotient by the Z 2 symmetry to leave just one copy of the 
desired space-time. 

In order to support the discontinuity at the IR wall and thus solve the equations of motion 
at the wall, there must be a source of stress-energy at r = r^. Therefore we will add an 
action, S wa u, localized to the wall and solve the equations of motion. One way that this 



is different from the situation discussed in 27 is that while those authors needed only to 



add a localized cosmological constant to the wall (as everything was only a function of the 
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radial variable), we now have spatial (x,y) dependence, so our boundary action must also 
have spatial dependence, S waa = S wa u(x,y). 



Which terms in the equations of motion will contribute to the boundary stress-energy? 
When integrating the equations of motion across the wall, the first derivative of any function 
of r will not contribute, whereas the second derivative will: 

drf(r) = /(r + e) - /(r - e) = 0; (3.30) 

ro—e 

drf"(r) = /'(r + e) - f'(r - e) = -2/'(r - e). (3.31) 

ro—e 

Therefore, in order to solve for the action at the wall, we only need to consider the terms 
in the equations of motion which have second derivatives of functions of r. At zeroth order 
in e the gauge field is independent of r, and the Einstein equations only depend on up to 
first derivatives of the metric functions. In this case, integrating the equations across the 
wall we find no contributions, and we find that we do not need an S wa u at zeroth order in 
e. 

At first order in e, we need to consider integrating the ip equation of motion across the 
wall. We will add the term 



SLu= J d*xV-h SmiW (3.32) 

r=r 

to the action, where h^ v is the induced metric at r = r , and 5m w is a localized shift in 
the mass of ip. The nonzero contributions to the equation of motion when integrated over 
the wall are 

/•ro+e pro+e 

— / dr^/^gg rr il)" = I dr\J—h Srn^ijjd^ — r ), (3.33) 

Jro-e Jro-e 
2 



which gives the result bm w = 2a/ L. 



Note that after adding the wall-localized mass term (3.32), the strategy for finding the 
critical field at which a phase transition occurs is the following. For a fixed choice of the 
wall localized mass and the location of the wall ro, there is a critical value of the 5-field 
at which the purely normalizable solution for tp obeys the boundary conditions. In this 
paper, we are always expanding about this critical field, with e parametrizing the distance 
from criticality. 

We note also that we will need to add additional terms to S wa ii when we consider the 
equations of motion at 0(e 2 ) in the next section. 
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3.4 Higher order corrections to the gauge field and metric 



The relevant equations of motion are the Einstein equations and the Euler-Lagrange equa- 
tions for ipjAfj,, equations (3.5) and (3.6). 

In an AdS2 x R 2 background, all the unknown functions scale as power laws in r. At 0(e) 
there is only the ip equation of motion. From §3 we know there exists a lattice solution of 
the form 



tpi(r, x, y) = p r a > e u i e 2 v n 

l=—oo 

where po is the magnitude of ipi and the scaling exponent is 

a = X - (l + v/l + 4(Q + m 2 L 2 )) 



(3.34) 



(3.35) 



■01 acts as a source in the gauge field equation of motion and Einstein equations at 0(e 2 ). 
Therefore we can extract the r scaling in the 0(e 2 ) corrections and solve the equations of 
motion for the spatial dependence. We write 

f t (r,x,y)=p 2 r 2a f t (x,y) (3.36) 

where fi = a, b, af , a|. By assuming a normalizable radial dependence of this form for each 
field, we are implicitly setting one integration constant to zero per function. Our choice of 
solution for each of these fields will also fix the form of the localized stress energy we will 
need to add at the wall in order to have a consistent solution. 

Now we examine the differential equations at 0(e 2 ). The equation of motion for A r gives 
us the constraint 

d x a x 2 + d y a\ = 0. (3.37) 

Besides this, we have 2 additional gauge field equations of motion (one each for x, y) and 5 
nontrivial Einstein equations (for G tt ,G rr ,G xx = G yy ,G rx ,G ry ) at 0(e 2 ). These are seven 
equations and four unknown functions. Luckily, three of them are redundant and we can 
find a consistent solution once we have chosen the form of the source, The Einstein 
equations are 

2(d y a x 2 -d x a y 2 ) + Q(d 2 x + d 2 )(a + b) + 2Q(4a 2 -1)6 = S^i) 
2(d x a y 2 -d y a x 2 ) + Q(d 2 x + d 2 )(a-b) + 2Q(2a + l)b = S 2 {^) 
2aa v 2 + Q(a — l)d x a — Qad x b = S^ipi) 
2aa 2 — Q(a — l)d y a + Qad y b = S^tjji) 
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2{d x a\ - d y afj - 2Q{a - l)(2a - I) a + 2Q(2a 2 - a + 1)6 = S 5 (^i) (3-38) 
and the gauge field equations are 

(d 2 x + d 2 + 2a(2a - l))a x 2 - Qd y b = S 6 (V>i) 

(8g + ^ + 2a(2a-l))a| + Q9 a> 6 = SV(V>i), (3.39) 
where the V'-dependent source terms are given by 

Si(i/>i) = "|(2a 2 + Q 2 m 2 + 2QV)|^| 2 + iQ 2 y(tx^i ~ V>i^D ~ Qd^iT + \d y ipi\ 2 ) 
S2W1) = § (~2a 2 + Q 2 m 2 + 2Q 2 y 2 )\^\ 2 - iQ^yWMi ~ 1>id*Pi) + Q(l^i| 2 + l^i| 2 ) 

SsM = ^\H 2 

Wi) = -^i^il 2 

SbM = "|(2« 2 + Q 2 m 2 + 2Q 2 y 2 )\^\ 2 + zQ 2 y(VWi - V>i^D - Q(|^i| 2 - \d y ^\ 2 ) 

Setyi) = -lQ 2 (rid x ^i-^id x ri+2zyQ\H 2 ) 

= l -Q 2 (i>id v ri-ridM (3.40) 

and a, b, af,a|,^i are now only functions of as we have omitted the power law r- 
dependence. 

We know that the vortex lattice solution is periodic in x, y with periodicity ^= in the x 
direction and v 2 ^q in the y direction (this is only for the rectangular lattice); therefore we 
can expand each of these functions as a double Fourier series in x, y as 

k,l 

where fi = a, b, af, a\, and we have pulled out the exponential function of m, n which will 
be present in all of the source terms. Notice that the periodicity implies that each unit 
cell has a net flux density of 7?. It remains to Fourier transform the source terms in the 



equations of motion in order to bring them into the form of equation (3.41 ), and then solve 
algebraic equations for the polynomial coefficients fi(k,l). In order to do this we will use 
properties of exponentials and the Fourier transform to write an infinite sum of Gaussians 
as an infinite sum of exponentials, 



k k v n 
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First we will do this for Q = 1. In this case we also have L = l/y/2, A = —1, and 



m 2 = 2(a 2 — a — 1). Plugging in our ansatz of equation (3.41), we get the following 
algebraic equations for the fi(k,l): 



2i ( ZvitTf a\ ) + 2(4ct 2 - 1)6 

' 2i\k ~ 

2/ I — a\ - lv x a\ ) + 2(2a + 1)6 



2-Kk 
2irk 



+ (Ivif 



{a + b) 
(a-b) 



2aa v 2 H — (a(a — 6) — a) 

2aa\ + ilv\(a — a (a — 6)) 



2? 



2nk ~„ 



a\ - lv x a% - 2(a - l)(2a - l)a + 2(2a 2 - a + 1)6 



2a(2a - 1) - 



2a(2a - 1) - 



2nk 

Vl 



2-nk 

Vl 



- (ivif 



{IVif 



af — ilvib 



al + 



2-nikr 
6 



a 



(l-2a) 7r 3 / 2 fc 2 _^|_ 

+ + 4y/7l 



2^ 



2v^ 

ia^/irk 
2v~i 
ilavi 

~ 

1 + a - 2a 2 

27^ 

ilvi 



3 / 2 k 2 l 2 v\ 



4A 
ik\fn 



2vi 



(3.43) 



From this we can see that we expect a and 6 to be real and af and 0% to be pure imaginary. 
For k = I = we get the solution 

a(2a 2 + a-4)-l ~ a 



" 4(a-l)(4a 2 -l) v /F " 4(2a + l)v^' °* ° 2 °' 
and in all other cases the solutions are 

av 2 {{a + l)4fc 2 7T 2 + v 2 {2 + {a + l)l 2 v 2 - 2a(2a 2 + a - 4))) 



(3.44) 



a(fc, /) = 
b(k,l) = 

af(M) 



al(k,l) 
where 



D 

16fc 4 vr 4 + (8k 2 7T 2 v 2 + 2l 2 vl)(l + l 2 v\ + a(2 - 3a)) + v 4 (-/S 4 + 4a 2 (2a 2 - 3a + 1)) 



2D 



ilv\{Ak 2 7T 2 + w 2 (l + Z 2 i; 2 + (2 - 3a)a)) 



2vrzA;t; 1 (4A; 2 7r 2 + w 2 (l + Z 2 v 2 + (2 - 3a)a)) 
D 



D = 2^(16k\ 4 +8k 2 7T 2 vf(l 2 vi+2a(l-2a))+vt(l%t+Al 2 vfa(l-2a) 

(3.46) 



(3.45) 



Note that the equations of (3.43) are not solvable for a = ±|, 1. In the case we have 
chosen, where (J = 1, a = |(1 + V5 + 2m 2 ). Thus, we can only solve these equations for 
some values of m. 
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3.5 0(e 2 ) stress-energy at the wall 



At 0(e 2 ) we need to consider the gauge field equations of motion and the Einstein equations 
integrated across the wall. We will now consider the following action at the wall: 

S wM = f d 3 xV=h{5m 2 M 2 + 6 2 A^ + 6 2 (T w )»} : (3.47) 

r=r 

where we have added a current which couples to the gauge field, as well as a source of 
stress-energy (T w )^ v localized at the wall. This is the most general form of action we can 
add to the wall and should easily lead to a solution. Because we don't want the boundary 
current or stress tensor to enter into the equations of motion at zeroth order, we have 
assumed that each term enters the action at 0(e 2 ). 

First we will consider integrating the gauge field equations of motion across the wall. The 
relevant equations are those for A x ,A y . The equations we must solve are 

/ drV=gg rr 9r(aT' v )" = / dr^h xx (J w ) x , y S(r-r ), (3.48) 
which have the solutions 

Act 

(Jw)x, y = -—a x ' y (r ,x,y). (3.49) 

Finally, we must consider the Einstein equations. There are three equations which include 
second derivatives of the fields, G tt , G xx , and G yy . The total stress-energy from the action 
at the wall takes the form 

V " T ' h ^C wall + e 2 ^h{A^{ J w ) u + (7^) (3.50) 



2 

where C w is the integrand of S wa u. After integrating the Einstein equations, we get the 
following set of equations for T w : 



a 



b(r ,x,y) = ^ (sm 2 w \Mro,x,y)\ 2 + jj(J w ) x + j^(T w ) tt + ^(T w ) xx + ^(T w ) yy ^j 

a L 5 ( Qy r 2 1 1 

— (a(r ,x,y) - b(r ,x,y)) = -— ( 5m 2 w |^| 2 - jj(J w ) x - j^(T w ) u - j^{T w ) xx + — 

^(a{r ,x,y) - b(r ,x,y)) = (sm 2 w |^i| 2 + jj(J w )x ~ j^(T w ) tt + ^(T w ) xx - ^ 

which have the solution 

(T w ) tt = -^(a(r ,x,y) - b(r ,x,y)) + ^^i(r ,x,y)| 2 
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Figure 1: The scalar vortex lattice configuration ipi(x,y). 



(T w ) X x — (T w ^)yy + a 2 (^0) •£? y) 

{T w ) yy = -2aL|^i(r ,2;,?/)| 2 + ^(36(r , x, y) - a(r Q ,x,y)). (3.51) 



We note that, as with the original Randall-Sun drum matching 27 , the wall-localized 
stress-energy violates the Null Energy Condition. This is not a significant concern here (as 
it was not there); warped solutions microscopically realizing Randall-Sundrum like warping 
have been found in the full string theory, and we expect similar solutions could be found in 
this more involved case. It does mean that the wall should not be considered as a 'brane' 
which has Goldstone modes that allow it to fluctuate in the transverse dimensions. 



3.6 Pictures of the modulated phase 

We conclude this section with representative plots of the scalar supporting the vortex 
lattice ipi^x^y) (Figure 1), the modulation of flux density in the crystal (Figure 2), and a 
representative crystalline metric function (plotted as a function of (x,y) in Figure 3 and 
(r,y) in Figure 4). All plots are for values of the parameters given by: Q — 1, a — | + 
Vi = a/2~7t (a square lattice). The functions have been approximated keeping 121 terms in 
the Fourier series (i.e., with k, I running from —5 to 5 in the formulae above). 
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Figure 2: The 0(e 2 ) correction to the A x gauge field, which controls the modulation of the 
flux density in the lattice. 




Figure 3: The metric function b(x, y) generated by the back-reaction of a square vortex 
lattice. 
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Figure 4: The metric function in front of dx 2 + dy 2 , now plotted as a function of y and r. 
We have chosen x = 2 for this plot. 



4 Comments on physics of the lattice model 



From the form of the deformed metric in §3, we can infer some basic facts about the 
physics of the lattice solution. The IR wall geometry we have implemented is a bottom-up 
implementation of IR confinement |28|. In physical observables, powers of the IR radial 
cutoff r can be replaced by powers of 1/A, with A the scale of confinement. However, it 
is common in such solutions that also at finite temperature, one could (after the transition 
from confinement to deconfinement represented by a horizon at some r < 7*0) replace 
powers of tq by 1/T. Using this correspondence, we can infer the leading corrections to 
thermodynamic functions. 

The free energy density J 7 will receive a correction at 0(e 2 ). It will have the general form 

7 ~T (l + e 2 T~ 2a + •■■) (4.1) 

where the leading term comes from the AdS2 geometry (and gives rise to the notorious 
extensive ground-state entropy), and the sub leading term is due to the physics of the vortex 
lattice. One can see that the 0(e 2 ) corrections will scale like A~ 2a in the confining geometry 



quite explicitly, both from the form of the wall action (3.47), and from the e expansion of 
the contributions to the bulk action. 



The schematic formula (4.1 ) makes it clear that for a given value of e and a, there is an IR 



scale beneath which one should not trust perturbation theory. To avoid this region, one 
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must keep 

T > e» . (4.2) 

As a increases, the regime of trustworthiness of the linearized solution shrinks; this is in 
keeping with the simple intuition that the perturbation expansion in powers of er a will 
break down at smaller values of r for larger a. 

Free classical defects would contribute a correction to the free energy density proportional 
to T (and, of course, inversely proportional to the lattice spacing). The exponent a there- 
fore parametrizes an anomalous scaling of the free energy per vortex, characteristic of the 
strongly coupled field theory. 

What happens beyond the regime where perturbation theory around the transition is valid? 
One natural speculation is that as one proceeds to the deep IR, the different lattice sites 



'decouple' in a manner similar to that seen in AdS% fragmentation 29 . It is possible 
that this would proceed via another phase transition (at a temperature/energy scale lower 
than the transition to the vortex lattice state) to a 'fragmented' state. Such a fate was 
proposed in |24] for the D-brane lattice models of [2], where it was speculated that this 
might also characterize the physics of generic AdS% horizons. The growing localization of 
the dominant contribution to the low-temperature entropy on distinct lattice sites in the 
gravity solution provides support for this idea, in perturbation theory. 

Finally, it is worth emphasizing that the lattices discussed here are quite distinct from those 
obtained in related literature by considering a periodic spatial variation in the chemical 
potential fi(x) [lO] . The key difference lies in the nature of the IR behavior. In systems 
with a finite charge density, spatial modulations of fi can and will be cancelled by the 
background charge carriers - they will be screened. The hard lattices of the sort discussed 
here, in contrast, cannot be screened (physically, one cannot screen a magnetic field), and 
their effects should be expected to persist to the deep IR. In the S-dual perspective, such 
a feature is natural for the analog of 'Wigner crystallization' of charged carries that are 
added to a conformal field theory. 



5 Connecting with more general gravity solutions 



Here, we describe how the lattice solutions we found in §3 should also arise in 'IR com- 
pletions' of metrics with rather general dynamical critical exponent z and hyperscaling 



violation parameter 9 30 33 . The basic point will be that, as in 34 and 35 , the AdS% 



can arise in the deep IR, where corrections to the action supporting such solutions can 
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become important. 



The bigger physics picture is the following. As discussed in 16 , one can expect expectation 
values of monopole operators to serve as order parameters for translation-breaking phases 
in doped critical field theories. By S-duality in the 4d bulk, one can map the magnetically 
charged field dual to the monopole operator, to an electrically charged field. The doping 
maps to a background magnetic field. Then, the lattices found in §3 give concrete examples 
of the solids described in [16] , in strongly coupled quantum field theory. The considerations 
of this section show that this can happen in models with rather general z and 9. 



5.1 Basic EMD theory and magnetic solutions 

We start with the bulk gravity theory represented by an Einstein-Maxwell-Dilaton action 
S = J d A x^g (R - 2(«90) 2 - V" - V (</>)) (5.1) 
where the gauge-coupling function is of the form 

f{4>) = ^ (5-2) 

and the scalar potential takes the form 

V&) = ^e<* . (5.3) 
This theory supports solutions of the form 

ds 2 = L 2 (-airfdt 2 + + b(r) 2 (dx 2 + dy 2 )] (5.4) 
V a ( r ) J 

with scalar profile 

(j)(r) = K log(r) . (5.5) 

In the simplest solutions, a and h take power-law scaling forms, and one finds 

a( r ) = C a r 1 " 7 , b(r) = C b r p (5.6) 
in the vicinity of the horizon (with differences arising as one goes towards the UV, if one 



wishes to find asymptotically AdS solutions, as in e.g. 30 ; note that we use the convention 
that r — > oo is the UV and r — > is the IR in this section, in contrast with earlier sections). 
The two exponents 7 and f3, fixed by the parameters a and 77 in the action, capture the 



20 



scaling properties of the IR fixed point induced by doping the CFT. The physically relevant 
dynamical critical exponent z and the hyperscaling violation exponent 9 are determined in 
terms of 7 and (3 by the equations 

2 



General values of these exponents were first obtained in dilatonic systems in 31 



In many cases that have been studied, these solutions can in fact be supported in two 
different ways. If one studies an electrically charged black brane, Gauss' law yields the 
solution 

F = Q^^dt A dr . (5.8) 

One then finds extremal solutions where K < and </> — > 00 near the horizon r ~ 0. This 
means that the coupling is vanishing. As discussed in [30] , in the very near- horizon regime, 
the solution is then unreliable; in a full UV complete theory like string theory, higher 
derivative corrections will usually become important, because new light states appear as 
g = e Q * — > 0. This difficulty can be avoided by turning on a small temperature, since this 
cuts off the running of the dilaton; and the near-horizon solutions for finite T are simple 
to write down as well. 

However, in a 4d bulk, one can also use bulk electric-magnetic duality to find a represen- 
tation of the solution in terms of a magnetically charged black brane, i.e. a field theory 
immersed in a background magnetic field. This allows us to make contact with our discus- 



sion in §2 and §3, and with the picture of 16 : 



Suppose one is interested in studying the physics of monopole operators to diagnose the 



phase structure of the 'electric' model. One could introduce monopoles into the theory (5.1 ) 
and compute their correlators using semi-classical techniques in a multi-soliton background. 
However, it is easier to realize that by electric/magnetic duality, one can represent the 
monopoles as quanta of fundamental electrically charged fields in a dual theory, where the 



electric background (5.8) is dualized to a background magnetic field. 



• As mentioned above, the running dilaton indicates an 'IR incompleteness' of the solution 
- as the dilaton runs to extreme values, new corrections typically become important and 
deform the solution. For magnetically charged black branes in these dilatonic system, one 
possible result of the corrections is the emergence, in the deep IR, of an AdS2 geometry. 



This was discussed for Lifshitz scaling metrics in 34 and for general 9 and z in 35 . 
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The end result is that in critical theories dual to dilatonic systems with fairly generic z and 
8, if we are concerned mostly with the physics at very low energies, we can study monopole 
operators by considering the dynamics of electrically charged scalars in an AdS2 throat 
supported by magnetic flux. This provides a rather general setting where our analysis in 
§3 could be of relevance. 



6 Discussion 



There are many interesting directions for future exploration of the analytical vortex lattice 
solutions described here. We briefly mention some of these now. 



• It should be possible to find analogous perturbative crystalline geometries emerging 
directly out of solutions with various values of z and 9, without invoking the transition to 



an AdS? space-time 36 



• It would be natural to explore replacing the IR Israel thin wall considered here, with a 
black brane horizon. 



• One would like to compute simple correlation functions in these backgrounds. For in- 
stance, quasi-universal features have been seen in the transport properties of simple holo- 



graphic lattice models in 101. Their analogues in this system are worth exploring 36 



• Most ambitiously, it would be nice to find the full non-linear solution to the coupled 
set of partial differential equations that characterize the system. This would most likely 
rely on powerful numerical techniques. This program should yield new insights on the 
'fragmentation' phenomenon, and the eventual emergence of a solid in a 'confined' phase 
of the boundary gauge theory. 
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